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Abstract 



We study the region of small transverse momenta in qq- and gg-initiated 
processes with no colored particle detected in the final state. We present the 
universal expression of the 0{cq) logarithmically-enhanced contributions up 
to next-to-next-to-leading logarithmic accuracy. From there we extract the 
coefficients that allow the resummation of the large logarithmic contributions. 
We find that the coefficient known in the literature as is process depen- 
dent, since it receives a hard contamination from the one loop correction to 
the leading order subprocess. We present the general result of B^ for both 
quark and gluon channels. In particular, in the case of Higgs production, this 
result will be relevant to improve the matching between resummed predictions 
and fixed order calculations. 
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The process in which a system of non strongly-interacting particles of large invariant 
mass Q 2 (lepton pairs, gauge boson(s), Higgs boson, and so forth) is produced in hadronic 
collisions is a well studied subject in perturbative QCD M. At transverse momenta q\ 
of order of Q 2 the cross section can be computed by using the standard QCD-improved 
parton model. When qx becomes small the simple perturbative picture is spoiled. This 
happens because large logarithmic corrections of the form log ^- arise due to a non complete 
cancellation of soft and collinear singularities between real and virtual contributions. These 
large logarithmic corrections can be resummed to all orders by using the Collins-Soper- 
Sterman (CSS) formalism 0. 

We consider the class of inclusive hard scattering processes 

Mi2-> A 1 + A 2 ...A n + X (1) 

where the collision of the hadrons hi and h 2 produces a system of non strongly-interacting 
final state particles A\ . . . A n carrying total momentum Q and total transverse momentum 
qr- According to the CSS formula, and neglecting terms which are finite in the limit qx — > 0, 
the cross-section can be written as [] : 

da ^ r 1 , r 1 , r.,L„ > da [ r LO) 



dq^dQ 2 



V I dxi f dx 2 f db-J (bq T ) b{Q 2 - x x x 2 s) 

~^Jo Jo Jo 2 d(p 

(f a/hl ® C ca ) L, |) (f b/h2 ® C- Cb ) L, |) S C {Q, b) , (2) 



where dcp = dPS(Q — > qi, q%, ■ ■ ■ q n ) represents the phase space of the system of non-colored 
particles, bo = 2e~ le and is the leading-order cross-section (i.e., with no final state 

partons and therefore qr = 0) for the given process (c, c can be either or g, g). The 

function C ab in Eq. (g) is a process-dependent coefficient function, Jo(bq T ) is the Bessel 
function of first kind and fi/h corresponds to the distribution of a parton i in a hadron h. 
The large logarithmic corrections are exponentiated in the Sudakov form factor 



$c(Q, b) = exp \ - / 

Jb, 



2 



dq 

byv q 2 



Q 



2 



A c (a s (q 2 )) \n^- + B c (a s (q 2 )) 
q 2 



(3) 



The functions A c , B c and C ab in Eqs. (fJH) have perturbative expansions in 



o:. s 



n=l 



M* s ) = E ( g ) 4 n) , (4) 



n=l 



BcM = E ( . (5) 



n=l 



C ab (a s ,z) = 5 ab 5(l - z) + E ( C£'(s) • (6) 



x It is assumed that all other dimensionful invariants are of the same order Q 2 . 
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In order to obtain the coefficients in Eqs.(|4]-|5|) at a given order, the differential cross- 
section at small qx has to be computed at the same order. A comparison with the power 
expansion in a s of the resummed result in Eq. (0) allows to extract the coefficients that 
control the resummation of the large logarithmic terms. 

In this letter we study the behaviour of cross-sections at small transverse momenta at 
second order in a s both in the quark and gluon channels. We find that the analytic form of 
the logarithmically-enhanced contributions can be computed perturbatively in a universal 
manner by using the recent knowledge on the infrared behaviour of tree-level and one- 
loop |5j QCD amplitudes. In this way, we are able to extract the coefficients A£\ Bjp, 
A^f> and B^ for any qq or gg initiated process in the class ([I]). Details on our calculation 
will be given elsewhere ||. Here we only present and discuss our main results. 

By following Ref. |7J we multiply the differential cross-section, calculated at parton level, 
by q T and take moments with respect to z — Q 2 / s defining the dimensionless quantity: 



/N QtQ 
dZzN da,/d<pdqWd<p- (?) 

In the quark channel, for the sake of simplicity and in order to compare our result for S(iV) 
to the one originally obtained for Drell-Yan in Ref. H], we restrict our attention to the 
non-singlet contribution to the cross-section defined by 

//' 

To have qx ^ at least one gluon must be emitted, thus E(iV) has the expansion: 

S (A0 = E!l £«(jv) + ( ^±)\m( N ) + . . . ( 9 ) 

2lT V27T/ 

In the following we will systematically neglect in £(iV) all contributions that vanish as 
q T -> 0. 

In order to compute the small qx behaviour of £(iV) our strategy is as follows. The 
singular behaviour at small qx is dictated by the infrared (soft and collinear) structure of 
the relevant QCD matrix elements. At 0(a s ) this structure has been known for long time 0. 
Recently, the universal functions that control the soft and collinear singularities of tree-level 
and one-loop QCD amplitudes at 0(a 2 s ) have been computed @,|]. 

By using this knowledge, and exploiting the simple kinematics of the leading order sub- 
process, we were able to construct improved factorization formulae that allow to control all 
infrared singular regions avoiding any problem of double counting 0. We have used these 
improved formulae to approximate the relevant matrix elements and compute the small qx 
behaviour of E(iV) in a completely universal manner. 

The calculation at 0(a s ) is straightforward and we recover the well-known results: 

sg(iv) = 2cv log % - 3cv + 2 7 ;; ) (ao (10) 

q T 

and 
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SW(iV) = 2C A \og % - 2(3, + 2 T W(AT) . (n) 
It 



Here C F = N £ N 1 , Ca = iV c and Tr = 1/2 are the SU(N C ) QCD colour factors, /?o = 
^C^ — and 7^(iV), J^y{N) are the quark and gluon one- loop anomalous dimensions, 

respectively. From Eqs. (TO,ffl) one obtains: 

A^ = 2C a B^ = -2 la a = q,g (12) 

where C a and 7 a are the coefficients of the leading (1 — z)~ l singularity and 5(1 — z) term 
in the one-loop Altarelli-Parisi kernels P aa , respectively, 

C q = C F C g = C A l q = \c F l g = Po- (13) 

At this order it is possible to obtain also the coefficient C^ b by considering the q F inte- 
grated distribution and including the renormalized virtual correction to the LO amplitude 
cc — > Ai + A 2 . . . A^, summed over spins and colours, which, at O(e ), can be written as0 

M2» MAIS M + cc. = £ (^fj (-^ - ^ + AW) |*S>(#)|'. (m) 

In Eq. ( |14"D the structure of the poles in e = (4 — d)/2 is universal || and fixed by the 
flavour of the incoming partons. The finite part A (which can depend on the kinematics of 
the final state non-colored particles) depends instead on the particular process in the class 
(ED) we want to consider. In the case of Drell-Yan we have 151 : 



A° Y = C F (-8 + ^) , (15) 



whereas for Higgs production in the mt op — > oo limit the finite contribution is |10 



Af = 5C A + lc A 7r 2 -3C F = 11 + 2n 2 . (16) 
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By using the information in Eq. (|TJ) we obtain for C t 



(i) 

ab 



cffto = -P e M + - z) (c a £ + \M4>)} (17) 

where Pab( z ) is the 0(e) term in the Altarelli-Parisi P a b(z, e) splitting kernel, given by: 



2 A11 our results are obtained using the factorization and renormalization prescriptions of the MS 
scheme and within the framework of conventional dimensional regularization. 
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-c F (1 - 

-C F z 

-2T R z{\ 

0. 



(18) 



At order a s the coefficients A)p and are fully determined by the universal Altarelli- 
Parisi splitting functions. The function C$ depends instead on the process through the one- 
loop corrections to the LO matrix element. The general expression in Eq. (|ITD reproduces 



correctly the coefficient Cj£ f° r Drell-Yan J7|, Higgs production in the m top — > oo limit fTl 



and 77 production [ffS 

At second order in a s , two different contributions to E( 2 )(iV) have to be considered: the 
real correction corresponding to the emission of one extra parton (i.e., two gluons or a qq 
pair) with respect to the 0(a s ) contribution, and its corresponding virtual correction. 

The double-real emission contribution is the most difficult to compute. One has to 
integrate over the phase space of the two unresolved final state partons keeping q T fixed and 
finally perform the z integration in Eq. (^). We find that, likewise YP-\N), this contribution 
to £( 2 )(iV) is process independent, i.e., it does not depend on the particular process in the 
class ([!]) we want to consider. 

The virtual contribution is simpler to compute and we find it to be process dependent. 
More importantly, its process dependence is fully determined by the function A appearing 
in the one-loop correction to the LO subprocess (see Eq. fll4])). 

In the following, for the sake of simplicity, we present the total results for YP>(N) 
corresponding to the choice of the factorization and renormalization scales fixed to Q 2 . 
Since we are interested in extracting the coefficients A^ g and -BPJ, as in the 0(a s ) case we 
concentrate on the diagonal qq and gg contributions to X( 2 )(iV). 

In the quark (non-singlet) channel we obtain: 



sg } (iv) 



lo s — 
+ 10 § — 



-2C 2 F 



9C 2 F + 2C F f3 



-IT 



7)+C F C A 



+ (2P + 12C F ) 1 ^(N)-A( 1 ^(N) 



-C F n f T R + 2C F A q {(f>) 



+ 4C 2 



+ 



cl 



C f Ca ( — 



13 11 



-7T 



+ 6C(3) 



+C F n f T R 



, 2 - 

1 + -7C' 
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+2 7( (2) ) (iv)+2c F 7;; ) (^) 



(N + l)(N + 2) 2 / 
- ZC F A q {(j>) 

1 \ 



+ 2 



{N + l)(N + 2) 



(i) 



The actual expression of the coefficient C qq 
correct. 



for ZZ production reported in Ref. [13] is not fully 
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+2j$(N)A q (<t>) - 2Cf(/5 + 3C F ) 



whereas for the gluon channel the result is: 



(N + l)(N + 2) 



r 



(19) 



99 



3Q 2 











8 


2 







C 



-2C\ 



+ \og 2 \ 8C A /3 -6C Al $(N) 



It 



'67 
9 + 3 



20 
~9~ ( 



+ 



C 2 



+2(3 ( 7 W(JV) - /%) - 4 ( 7 £>(JV) - ft) 2 - 4n /7 «(iV) 7 «(iV) 
A I -y + 2C(3)) + 2C F n f T R + *-C A n f T R - 2(3 U^) + C A ^\ 



7T 



+2y%(N) + 2j£(N) [AM + Ca~)+ ±C F n n ^(N)-^—^ 



(20) 



In Eq. (0) 7$(iV) i s the non singlet space- like two- loop anomalous dimension |]Tq| , in 
Eq. ( pD| ) 7^(iV) is the singlet space-like two- loop anomalous dimension |16| , is the 
Riemann ( function (C(3) = 1.202...) and the function A a (4>) is defined through Eq. (jHD . 
The coefficients —t , , on and (Ar i 2 i nave origin on the N moments of —P^Jz) and —P*(z), 



(N+l)(N+2) (AT+2) 

respectively. 

The iV-dependent part of the results in Eqs. ( |T9"l , pO"D agreesQ with the one obtained from 
the second order expansion of Eq. (fj) (see e.g. Ref. By comparing also the N- 

independent part we obtain for 



A® = KA& 



K = C, 



'67 

18 6 



n f Tr- 



io 
9" 



(21) 



in agreement with the results of Ref. |l7|Jl8| . Moreover we find that can be expressed 
as: 

= -2 8P$2 + Po ("O^ 2 + A(0)) a = q,g (22) 

where 8Pf£} are the coefficients of the 8(1— z) term in the two- loop splitting functions Pf£}(z) 



r5|jr6fl , and are given by 

3 7T 2 



J + 6C(3)) + C F C A (I + ^ - 3C(3)) - M (I + ^ 



<5P^ = C 2 A [- + 3C(3)J - C^T* - -C A n f T R . 



(23) 



^This can also be regarded as an independent re-evaluation of the two-loop anomalous dimensions. 
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From Eq. ( p2|) we see that B^ 2 \ besides the —2 5P^ term which matches the expectation 
from the 0(a s ) result, receives a process- dependent contribution controlled by the one-loop 
correction to the LO amplitude (see Eq. (0)). We conclude that the Sudakov form factor 
in Eq. (^) is actually process dependent beyond next-to-leading logarithmic accuracy. The 
interpretation of this result will be given elsewhere ||19|| . 

However, by using the general expression in Eq. ( p2"D it is possible to obtain B^ for 
a given process just by computing the one- loop correction to the LO amplitude for that 
process. For the case of Drell-Yan, by using Eq. (|i~5l), our result for 'Efq(N) agrees with the 
one of Ref . J7) , and we confirm: 

3 \ „„/!!, 193 „ w „,\ „ m /17 4 2 



Bf )DY = C 2 (V - ^ - 12C(3)) + C F C A (^ 2 - ^ + 6C(3)) + C F n f T R 



7T 

3 9 



(24) 



In the case of Higgs production in the m top — > oo limit, by using Eq. ([16]) we find: 



Bf )H = C\ (| + f vr 2 - 6 C (3)) + AC F n f T R - C A n f T R (| + - ^C F C A . (25) 

In particular, this result allows to improve the present accuracy of the matching between 
resummed predictions |20| and fixed order calculations [f2~T|j . 



Summarizing, we have studied the logarithmically-enhanced contributions at small trans- 
verse momentum in hadronic collisions at second order in perturbative QCD. The calculation 
was performed in an process-independent manner, allowing us to show that the Sudakov form 
factor is actually process dependent beyond next-to-leading logarithmic accuracy. We have 
provided a general expression for the coefficient for both quark and gluon initiated 
processes. 

We would like to thank Stefano Catani and Zoltan Kunszt for many valuable discussions 
and comments, Christine Davies for providing us with a copy of her Ph.D. thesis and Luca 
Trentadue, James Stirling and Werner Vogelsang for discussions. 
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